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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (met¬ 
ric)  units  as  follows: 


Multiple 


feet  per  second 

inches 

kips  (force)  per 
square  inch 

miles  (US  statute) 

pounds  (force) 

pounds  (force)  per 
cubic  foot 

pounds  (force)  per 
foot 

pounds  (force)  per 
inch 

pounds  (force)  per 
square  inch 


0.3048 

0.3048 

25.4 

6.894757 

1.609 

4.448222 

157.08748 

14.593904 

175.1268 

0.006894757 


To  Obtain 
metres 

metres  per  second 

millimetres 

megapascals 


kilometres 

newtons 

newtons  per  cubic 
metre 

newtons  per  metre 


newtons  per  metre 


megapascals 


SUMMARY  OF  THE  RICHARD  B.  RUSSELL 


CONCRETE  DAM  VIBRATION  STUDY 


PART  I:  INTRODUCTION 


1.  The  Richard  B.  Russell  Dam,  recently  built  by  the  US  Army  Corps  of 
Engineers  (CE),  is  approximately  170  miles*  from  the  mouth  of  the  Savannah 
River  between  Georgia  and  South  Carolina.  As  shown  in  Figure  1,  the  crest  of 
the  concrete  gravity  part  of  the  dam  is  1,884  ft  long.  It  is  composed  of  13 
nonoverflow,  8  intake,  and  11  spillway  monoliths,  the  tallest  of  which  is 
approximately  200  ft  high. 

2.  Results  of  prototype  vibration  tests  of  the  Richard  B.  Russell 
Dam,  with  and  without  reservoir,  were  compared  with  linear  elastic  three- 
dimensional  dynamic  finite  element  (FE)  analyses  of  the  dam  (Chiarito  and 
Mlakar  1983;  Bevins,  Chiarito,  and  Hall,  in  preparation).  This  comparison 
provides  a  basis  for  assessing  the  accuracy  of  linear  elastic  FE  models  for 
predicting  dynamic  properties  of  concrete  gravity  dams. 

3.  A  modal  assurance  criterion**  (MAC)  was  used  to  compare  the  results 
of  experimental  and  analytical  structural  dynamics  (Allemang  and  Brown  1982). 
(See  Appendix  B  for  a  discussion  of  the  modal  assurance  criterion.) 

4.  The  grid  for  the  three-dimensional  FE  analyses  of  the  dam  is  shown 
in  Figure  2.  Four  FE  analyses  were  conducted: 

a.  Fixed  base  or  foundation,  without  reservoir. 

b.  Fixed  base,  with  reservoir. 

c.  Flexible  base,  without  reservoir. 


d.  Flexiole  base,  with  reservoir. 


The  flexible  base  was  modeled  with  vertical  and  horizontal,  springs  using  a 
spring  constant  computed  from  material  properties  of  the  foundation.  Concen- 


:rated  masses  were  added  to  model  the  mass  due  to  the  reservoir  and  the 


tainter  gates. 


*  A  table  of  factors  for  converting  nor. -SI  units  of  measurement  to  SI 
(metric)  units  is  presented  or  page  3- 


Certam  technical  terms  are  defines  in  the  Glossary,  Appendix  E. 


PART  II:  PROTOTYPE  TESTS 


5.  Two  forced  vibration  tests  were  made  on  the  Richard  B.  Russell  Dam 
before  and  after  impoi  ldment  of  the  reservoir  to  determine  the  natural  fre¬ 
quencies,  mode  shapes,  modal  damping  ratios,  and  relative  joint  movements 
(Chiarito  and  Mlakar  1983;  Bevins,  Chiarito,  and  Hall,  in  preparation).  The 
first  low-level  forced  vibration  test  (test  1),  before  reservoir  impoundment, 
was  conducted  during  January  and  February  1982  with  a  mean  headwater  elevation 
of  393.  The  second  test  (test  2),  with  reservoir  impoundment,  was  conducted 
during  June  and  July  1989  with  a  mean  headwater  elevation  of  970.  These  tests 
provided  an  experimental  measure  of  the  prototype  hydrodynamic  interaction  and 
a  comparison  of  the  changes  of  the  dynamic  properties  produced  by  the 
reservoir . 

6.  The  dam  was  excited  at  monoliths  7,  16,  and  22  by  a  crest-mounted 
inertial  mass.  The  force  input  was  in  an  upstream-downstream  direction.  As 
shown  in  Figures  3  and  9,  respectively,  three  arrays  of  servo  accelerometers 
measured  the  horizontal  crest  accelerations  of  31  monoliths  in  the  first  test 
and  all  32  monoliths  in  the  second  test.  Measurements  from  the  drive  point  on 
monolith  16  overlapped  measurements  from  drive  points  on  monoliths  7  and  22. 

7.  In  each  test,  accelerometers  were  placed  at  different  elevations  in 
the  three  drive  point  monoliths  to  measure  the  horizontal  motions  in  the 
vertical  planes.  Figures  5  and  6  show  the  cross-sectional  accelerometer  loca¬ 
tions  for  the  first  and  second  tests,  respectively.  Pressure  gages  were  also 
placed  on  monolith  7  in  the  second  test  (Figure  6). 

8.  Relative  joint  motion  arrays  were  positioned  at  three  locations  in 
test  1  and  two  locations  in  test  2.  They  consisted  of  two  accelerometers 
closely  spaced  on  each  side  of  the  joint  being  monitored.  Tests  1  and  2  joint 
accelerometer  locations  are  shown  in  Figure  7. 

9.  The  abutments  of  the  sides  of  the  dam  restrained  the  end  monoliths 
with  some  degree  of  fixity.  However,  the  abutment  on  the  Georgia  side  was  not 
complete  during  the  first  test  and  was  approximately  60  ft  below  the  completed 
crest  elevation  A95.  Therefore,  a  portion  of  monolith  1  was  unrestrained  on 
che  Georgia  side  during  the  first  test. 

10.  Static  and  dynamic  measurements  on  6-  by  12-in.  concrete  cylinders 
from  the  dam  are  summarized  in  Tables  1  through  3-  Dynamic  modulus  of  elas¬ 
ticity  of  seven  rock  core  samples  are  listed  in  Table  9.  The  in  situ  elastic 


modulus  of  the  foundation,  ,  was  estimated  from  the  dynamic  modulus  of 

elasticity,  cross-hole  seismic  velocity,  and  pulse  velocity  of  the  rock  core 

samples.  As  discussed  in  Appendix  A,  the  estimated  value  for  E~  was 
£  * 

6.05  <  10  psi. 

11.  The  mode  shapes  measured  of  the  crest  of  the  dam  from  the  two  vi¬ 
bration  tests  are  shown  in  Figures  8  through  12.  The  response  at  monolith  22, 
a  drive  point  measurement,  was  neglected  because  it  was  an  order  of  magnitude 
larger  than  the  other  responses.  This  extremely  large  response  seemed  to  be 
in  error  and  distorted  the  mode  shapes.  By  visual  inspection  of  these  modes, 
the  responses  in  the  midsection  of  the  dam  compared  well,  but  differed  at  the 
ends  of  the  dam.  Slight  shifts  in  the  maximum  response  points  towards  the 
ends  of  the  dam  are  an  indication  that  the  hydrodynamic  effects  added  mass  to 
the  dam's  midsection.  Note  that  the  fourth  crest  flexural  mode,  Figure  11, 
was  determined  for  test  2  only.  This  is  because  the  vibrator  locations  may 
have  been  nodes  of  the  fourth  mode  in  test  1  (Chiarito  and  Mlakar  1983).  (The 
locations  of  the  nodes  were  possibly  different  for  tests  1  and  2.) 

’2.  As  expected,  the  reservoir  impoundment  decreased  the  natural 
frequencies  of  the  dam  due  to  the  added  mass.  The  damping  estimates  and 
naturae  frequencies  from  test  1  and  2  are  shown  in  Tables  5  and  6,  respec¬ 
tively.  With  the  exception  of  mode  5,  the  reservoir  had  little  effect  on 
the  damping  estimates.  Percent  decreases  in  natural  frequencies  were  the 
largest  for  mode  1  (10.1  percent),  with  smaller  decreases  for  the  higher 
modes  (Bevins,  Chiarito,  and  Hall,  in  preparation).  (As  discussed  later, 
natural  frequencies  from  finite  element  analyses  are  also  included  in  Table  6 
for  comparison.; 

'3.  The  MAC  was  used  to  compare  the  crest  mode  shapes  (Allemar.g  and 
Brown  1982).  (Mor.olith  22  was  not  used  in  computing  the  MAC  because  excessive 
local  responses  were  recorded  in  test  2  at  this  monolith.)  Table  7  summarizes 


the  results. 


Mode  3  showed  reasonable  agreement,  but  the  other  modes  had  fair 


to  poor  correlation.  Possible  explanations  of  the  poorer  correlation  between 
mode  shapes  include  the  following: 

a.  The  force  input  was  r.ot  iarge  enough  to  completely  excite  the 
entire  structure  an  the  lower  frequencies.  (A  typical  force  is 
shewn  in  Figure  '3).  Ideally,  two  or  more  force  inputs  would 
be  better.  By  visual  inspection  of  modes  1  and  2,  the  mid¬ 
section  of  the  dam  agrees  we .  ,  but  at  the  ends  the  shapes 
1  i  f •■‘►'•r  .  This  would  accent. t  :  the  low  MAC. 


b.  A  high  degree  of  modal  coupling  (closely  spaced  resonance 
frequencies)  existed  at  the  higher  frequencies. 

c.  The  friction  forces  between  monoliths  may  have  been  different 
for  test  1  and  2.  Tests  1  was  conducted  in  winter,  and  test  2 
was  conducted  in  summer.  Thus,  expansion  and  contraction  of 
the  material  could  have  changed  the  normal  and  frictional 
forces  between  monoliths. 

14.  Except  at  elevation  430  (normalized  height  was  0.65),  the  experi¬ 
mentally  measured  pressures  on  the  upstream  face  of  monolith  7  in  the  second 
test  with  reservoir  (Bevins,  Chiarito,  and  Hall,  in  preparation)  were  close  to 
theoretical  values  calculated  by  a  simplified  procedure  developed  by  Fenves 
and  Chopra  (’986).  (Due  to  a  scaling  error  in  the  gage  calibration,  an  ex¬ 
cessively  large  pressure  was  recorded  at  elevation  430.)  The  theoretical 
pressures  are  based  on  the  fundamental  cantilever  mode  shape  and  maximum  ac¬ 
celeration  of  a  two-dimensional  section.  The  measured  pressure  distribution 
also  was  in  agreement  with  pressures  measured  during  vibration  testing  of  the 
Pine  Flat  Dam  (Rea,  Liaw,  and  Chopra  1972).  The  magnitude  of  the  experiment¬ 
ally  and  theoretically  determined  hydrodynamic  pressure  for  the  fundamental 
frequency  (5.3  Hz)  of  the  dam  is  plotted  in  Figure  14  as  a  function  of  ele¬ 
vation  along  the  centerline  of  monolith  7.  (The  extremely  large  pressure  at 
elevation  430  is  not  shown  in  Figure  14.)  As  shown  in  Figure  14,  theoretical 
pressures  were  computed  for  two  concrete  moduli  to  establish  a  band  m  which 
the  experimental  pressures  should  fall.  This  is  because  the  monolith  was 
constructed  using  different  concrete  strengths  for  the  exterior  and  interior 
mass  (Tables  1  and  2).  Evidence  of  low  signal-to-noise  ratio  can  be  seen  in 
the  pressure  distribution  in  Figure  14  (i.e.,  the  points  do  not  form  a  smooth 
distr ibution ) . 

'5.  Measurements  from  both  tests  at  the  relative  joint  motion  arrays 
indicated  relative  motions  between  monoliths.  The  nonlinear  behavior  of  the 
dam  is  attributed  mostly  to  the  monoliths'  joint  behavior. 


PART  III:  FINITE  ELEMENT  MODEL 


'inite  Element  Program 


It.  A  structural  analysis  program  for  the  static  and  dynamic  response 
of  linear  systems  was  used  to  calculate  the  natural  frequencies  and  mode 
shapes  of  the  dam  (Bathe,  Wilson,  and  Petersor.  1973).  This  program  is  coded 
in  FORTRAN  66  and  was  run  op.  the  Honeywell  DPS  8  mainframe  computer  at  the  US 
Army  Engineer  Waterways  Experiment  Station: 

17.  The  program  contains  the  following  element  types: 

a.  Three-dimensional  truss  element. 

b.  Three-dimensional  beam  element. 

c.  Plane  stress  and  plane  strain  element. 

d.  Two-dimensional  axisynmetric  solid. 

e.  Three-dimensional  solid. 

f.  Variable-number-nodes  thick  shell  and  three-dimensional 
elements . 

2.  Thin  plate  or  thin  shell  element. 

h.  Boundary  element. 

i.  Pipe  element  (tangent  and  bend). 

These  structural  elements  can  be  used  in  static  or  dynamic  analyses  to  model 
a  large  number  of  two-  or  three-dimensional  structural  problems.  The  capacity 
of  the  program  depends  mainly  on  the  total  number  of  nodal  points  in  the  sys¬ 
tem,  the  number  of  eigenvalues  needed  in  the  dynamic  analysis,  and  the  com¬ 
puter  used.  Practically  no  restriction  is  placed  or,  the  number  of  elements 
used,  the  number  of  load  cases,  or  the  order  and  bandwidth  of  the  stiffness 
matrix.  Each  nodal  point  in  the  system  can  have  from  zero  to  six  displacement 
degrees  of  freedom.  Because  the  element  stiffness  and  mass  matrices  are 
assembled  in  condensed  form,  the  program  is  equally  efficient  in  the  analysis 
of  one-,  two-,  or  three-dimensional  systems. 

'8.  The  formation  of  the  str  tur°  m trices  is  carried  out  in  the  same 
manner  in  a  static  or  dynamic  analysis,  rear  typers  of  dynamic  analysis  can  be 
o r forme  i  by  t h n  p rograrr. : 

a.  Doternin.at  im-  •' f  y.:  a  men*1  shapes  and  frequencies  only. 
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d.  Dymamic  response  analysis  for  arbitrary  time-dependent  loads 
using  mode  superposition  or  step-by-step  direct  integration  of 
the  equations  of  motion. 

19-  The  structural  mode  shapes  and  frequencies  can  be  obtained  by  two 
distinct  algorithms: 

a.  The  determinant  search  method  is  used  for  small  problems  in 

which  the  total  stiffness  matrix  can  be  kept  in  the  core  of  the 

computer . 

b.  The  subspace  iteratio-’  method  is  used  for  large  problems  in 

which  only  portions  of  the  stiffness  matrix  can  be  retained  in 

the  core  of  the  computer. 


A  subspace  iteration  method  was  used  to  solve  for  the  undamped  natural  fre¬ 
quencies  and  mode  shapes  presented  in  this  report. 


Parameter  Studies 


20.  The  g^id  for  the  FE  model  was  developed  by  conducting  parameter 
studies  of  isolated  nonoverflow,  spillway,  and  intake  monoliths.  Modal  anal¬ 
yses  were  used  to  obtain  appropriate  FE  grid  sizes  for  each  of  the  monoliths. 

A  compatible  grid  was  found  for  all  of  the  monoliths,  thus  making  possible  the 
generation  of  the  nodal  and  element  data  with  the  computer  (Chiarito  1985). 
Favorable  results  were  obtained  when  modal  analyses  of  each  monolith  using 
uniform  grids  were  compared  with  analyses  using  the  grid  for  the  entire  three- 
dimensional  FE  model.  Figures  15  through  17  show  typical  sections  of  the 
monoliths  and  the  resulting  three-dimensional  grids  of  each  monolith.  (The 
resulting  three-dimensional  FE  grid  of  the  dam  was  shown  previously  in 
Figure  2 . ) 

21.  Figure  15  shows  the  penstock  tube  opening  in  the  intake  monolith 
modeled  as  solid,  or  "smeared-hole, "  elements  with  29  percent  less  stiffness 
and  density  than  the  mass  concrete.  Modal  analyses  showed  good  agreement 
between  a  finer  FE  grid  that  included  the  penstock  opening  and  the  grid  shown 
in  Figure  5  using  solid  elements. 

22.  Studies  were  also  performed  to  obtain  a  suitable  Winkler  spring 
representation  of  the  foundation.  Comparisons  of  two-dimensional  and  three- 
dimensional  FE  dynamic  analyses  on  the  nonoverflow  section  using  Winkler 
springs  with  analyses  using  an  FE  foundation  grid  and  a  fixed  foundation  are 
presented  in  Appendix  A. 
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Material  Properties 

23.  The  concrete  properties  used  in  the  FE  analyses  were  taken  from 
previous  measurements  on  6-  by  12-in.  concrete  cylinders,  as  summarized  in 
Tables  1  through  3.  The  concrete  samples  tested  in  1982  and  1989  were  5  and 
17,  respectively.  Generally  higher  concrete  strengths  were  measured  in  1982 
than  in  1984;  however,  one  cannot  conclude  that  the  concrete  was  stronger  in 
1982  than  it  was  in  1984  because,  statistically,  not  enough  samples  were 
tested  in  1982  to  indicate  any  definite  change  in  concrete  strength  between 
•982  and  1984.  For  example,  the  cylinder  strengths  for  the  3-ksi  design 
strength  concrete  determined  from  the  2  samples  tested  in  1982  fell  within  one 
standard  deviation  of  the  14  samples  tested  in  1984.  This  statistical  infor¬ 
mation  indicates  that  if  14  concrete  samples  had  been  tested  in  1982,  the 
range  of  data  would  have  been  similar  to  the  1984  results. 

24.  Dynamic  modulus  of  elasticity  values  used  were  5.1  *  10^  psi  for 
exterior  mass  concrete  and  4.26  *  10^  psi  for  interior  mass  concrete.  From 
previous  foundation  tests,  the  dynamic  elastic  modulus  of  the  foundation  used 
to  compute  the  spring  constant  was  6.05  *  10^  psi,  as  discussed  in  Appendix  A. 
Using  elastic  theory,  the  average  horizontal  and  vertical  spring  constant  was 
4  .  10 10  Lb /ft  in  the  flexible  base  model  (Timoshenko  and  Goodier  197 0). 

Element  Types 

25.  Four  element  types  were  used  to  model  the  dam.  Variable-number- 
nodes  18  to  21  nodes)  three-dimensional  solid  elements  were  used  to  model  the 
mass  concrete  of  the  nonoverflow,  intake,  and  spillway  monoliths  (Bathe, 
Wilson,  and  Peterson  1973).  Plate  eLements  were  used  to  model  the  concrete 
spillway  piers,  the  concrete  training  walls  at  spillway  monoliths  16  and  26, 
and  the  concrete  deck  of  the  spillway  bridge.  Beam  elements  were  used  to 
model  the  steel  girders  of  the  spillway  bridge,  and  Winkler  springs  were  used 
along  the  base  of  the  model  in  the  flexible  base  analysis. 

Boundary  Conditions 

26.  By  definition,  only  the  X  ,  Y  ,  and  Z  translations  were  defined 
on  the  three-dimensional  thick-shell  elements.  All  degrees  of  freedom  were 
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deleted  along  the  base  of  the  model  for  the  fixed  base  analysis.  For  the 
flexible  base  analysis,  the  longitudinal  (Y)  translations  and  all  rotations 
were  fixed  along  the  base  of  the  model.  Except  at  the  dam  abutment  interface, 
horizontal  and  vertical  springs  were  placed  at  all  nodes  on  the  base  in  the 


X  and  Z  directions. 


27.  In  order  to  model  the  restraint  due  to  the  abutments,  nodes  along 
the  two  end  sections  (monoliths  1  and  32)  had  all  degrees  of  freedom  deleted 
(i.e.,  nodes  fixed).  However,  a  few  of  the  nodes  had  translational  degrees 
of  freedom  defined  on  the  Georgia  side  for  the  model  of  the  dam  without 
reservoir.  As  was  discussed  previously,  this  abutment  was  not  complete  at  the 


time  of  the  first  vibration  test. 


28.  For  the  model  without  reservoir,  5,127  and  5,609  degrees  of  freedom 
were  used  for  the  fixed  base  and  flexible  base  FE  models,  respectively.  For 
the  model  with  reservoir,  5,100  and  5,582  degrees  of  freedom  were  used  for  the 
fixed  base  and  flexible  base  FE  models,  respectively. 


Modeling  of  Added  Mass 


29.  Concentrated  masses  were  added  at  nodes  thought  to  be  in  contact 
with  the  reservoir  to  model  the  reservoir  mass  and  the  tainter  gates.  A  con¬ 
stant  49.6-ft-wide  mass  of  reservoir  was  applied  to  the  model  with  reservoir. 
Studies  (see  Appendix  C)  showed  that  this  constant  width  of  reservoir  added 
the  same  amount  of  mass  as  a  modified  Westergaard  (1933)  analysis.  A  typical 
nonoverflow  monolith  for  the  dam  was  analyzed  with  a  constant  49.6-ft-wide 
mass  of  reservoir.  As  discussed  in  Appendix  C,  the  fundamental  frequency 
compared  well  with  the  result  obtained  from  an  approximate  procedure  by  Chopra 
(1978).  The  mass  of  the  reservoir  acting  on  the  tainter  gates  and  the  mass  of 
the  gate  was  applied  at  the  trunnion  locations  on  the  spillway  piers. 
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PART  IV:  FINITE  ELEMENT  MODEL  RESULTS 


Dam  Without  Reservoir 

30.  The  three-dimensional  shapes  (eigenvectors)  shown  in  Figures  18  and 
19  include  the  first  nine  major  bending  modes  of  the  dam  crest  computed  from 
the  FE  model  with  fixed  base.  At  7.7  Hz  (eigenvector  2),  the  major  response 
is  associated  with  the  sidesway  of  the  roadway  and  concrete  spillway  pier  sys¬ 
tem.  The  other  eigenvectors  show  that  the  major  responses  of  the  dam  are  as¬ 
sociated  with  bending  of  the  dam  perpendicular  to  the  longitudinal  axis  of  the 
dam.  The  fifth  eigenvector,  at  10.4  Hz  (i.e.,  the  fourth  bending  mode  of  the 
crest),  is  the  shape  computed  by  the  FE  analysis  but  not  measured  experimen¬ 
tally.  As  mentioned  previously,  it  is  possible  that  the  vibration  exciters 
were  located  at  nodes  of  that  mode  shape  during  the  first  vibration  test 
(Chiarito  and  Mlakar  1983). 

31.  Previously,  initial  three-dimensional  FE  analyses  were  conducted 
for  the  entire  dam  without  reservoir  (Chiarito  1985).  The  main  differences 
between  the  initial  analyses  and  the  latest  analyses  were  that  the  initial 
analyses  did  not  model  the  training  walls  at  monoliths  16  and  26  and  had  more 
defined  degrees  of  freedom  at  the  dam  abutment  interfaces.  Material  proper¬ 
ties  were  slightly  different,  and  the  initial  analyses  did  not  add  the  weight 
of  the  tainter  gates.  The  initial  analyses  had  natural  frequencies  approxi¬ 
mately  5  percent  greater  than  the  latest  analyses. 

32.  Also,  three-dimensional  dynamic  FE  analyses  were  conducted  on 
isolated  monoliths,  and  two-dimensional  FE  analyses  were  conducted  on  the 
nonoverflow  monolith.  Foundation  and  reservoir  effects  were  not  included  in 
these  analyses,  and  dynamic  modulus  of  elasticity  values  were  used.  The  fun¬ 
damental  frequencies  calculated  for  the  isolated  nonoverflow,  intake,  and 
spillway  monoliths  were  9.8  Hz,  6.9  Hz,  and  11.3  Hz,  respectively.  These 
results  indicate  that  the  intake  sections  are  the  most  flexible  part  of  the 
dam  and  that  the  spillway  monoliths  are  the  stiffest  part  of  the  dam.  Also, 
the  fundamental  frequency  for  the  two-dimensional  analysis  of  the  nonoverflow 
section  was  8.4  Hz.  Previously,  a  two-dimensional  analysis  of  the  nonoverflow 
section  using  a  static  modulus  of  elasticity  estimated  a  fundamental  frequency 
of  7.63  Hz  (Norman  1979) . 

33-  The  calculated  fundamental  frequency  of  the  entire  dam  without 
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reservoir  was  7.1  Hz  for  a  fixed  foundation.  This  frequency  was  less  than  the 
fundamental  frequencies  of  the  nonoverflow  and  spillway  monoliths  but  slightly 
more  than  the  fundamental  frequency  of  the  intake  monolith.  As  expected,  the 
intake  sections  increased  the  flexibility  of  the  entire  dam. 

Dam  With  Reservoir 

34.  The  three-dimensional  shapes  (eigenvectors)  shown  in  Figures  20 
and  21  are  the  first  eight  major  bending  modes  of  the  dam  crest  and  two  mode 
shapes  associated  with  sidesway  of  the  roadway  and  concrete  spillway  pier 
system.  These  eigenvectors  were  calculated  from  the  FE  model  with  fixed 
base.  The  eigenvectors  not  shown  are  also  associated  with  sidesway  of  the 
spillway  pier  system. 

Comparisons  of  Dam  With  and 
Without  Reservoir 


35.  Tables  8  and  9  summarize  results  of  the  FE  analyses  with  and  with¬ 
out  reservoir  for  fixed  and  flexible  bases,  respectively.  The  FE  mode  shapes 
with  and  without  reservoir  were  compared  using  MAC.  Most  mode  shapes  compared 
well,  except  at  the  higher  frequencies.  Figures  22  through  31  show  the  trend 
toward  poorer  correlation  among  FE  fixed  and  flexible  base,  mode  shapes,  with 
and  without  reservoir,  as  the  frequency  increases. 

36.  As  expected,  the  frequencies  computed  from  the  model  of  the  dam 
with  reservoir  were  lower  than  from  the  model  without  reservoir.  Table  6, 
presented  earlier,  is  a  summary  of  the  natural  frequencies  determined  by 
experiment  and  analysis. 


PART  V:  COMPARISON  OF  FINITE  ELEMENT 
MODEL  TO  PROTOTYPE  TEST  RESULTS 


Dam  Without  Reservoir 


37.  Comparisons  of  the  dam  crest  mode  shapes  are  shown  in  Figures  32 
through  36.  Note  that  the  mode  numbers  correspond  to  bending  modes  of  the  dam 
crest.  As  mentioned  previously,  the  fourth  mode  was  not  measured  experi¬ 
mentally,  but  was  determined  analytically.  Table  10  summarizes  the  results  of 
these  analytical  and  experimental  mode  shape  comparisons.  Based  upon  the 
results  of  the  MAC  values,  the  correlation  between  analytical  and  experimental 
mode  shapes  was  favorable.  Overall,  the  model  with  a  fixed  base  had  slightly 
better  mode  shape  correlation  than  did  the  model  with  a  flexible  base,  pos¬ 
sibly  because  the  base  of  the  dam  was  more  rigid  than  modeled  by  the  FE 
flexible  base  analysis.  For  a  qualitative  comcarison,  the  mode  shapes  all 
appear  to  have  the  largest  amplitude  in  the  same  general  location  on  the  crest 
of  the  dam.  The  MAC  does  indicate  quantitatively  hew  good  the  comparison  is 
among  modes. 

38.  As  shown  in  Table  6,  the  frequencies  estimated  by  the  FE  models 
were  higher  than  frequencies  measured  experimentally  ( 16  to  30  percent  higher 
for  the  FE  fixed  base  model,  and  9  to  18  percent  higher  for  the  flexible  base 
model).  This  is  reasonable  since  the  FE  models  are  approximations  of  the  dam 
using  a  discrete  number  of  degrees  of  freedom.  The  FE  model  in  this  case  is 
stiffer  than  the  actual  dam.  The  flexible  base  model  estimated  frequencies 
closer  to  the  experimentally  measured  frequencies  than  did  the  fixed  base 
model  because  foundation  interaction  decreases  the  frequency  (Norman  and  Stone 
1979). 

39.  Comparisons  of  the  experimental  mode  shapes  in  cross  section 
with  results  from  the  FE  model  with  fixed  foundation  are  shown  in  Figures  37 
through  56.  Note  that  the  frequencies  of  the  cross-sectional  modes  correspond 
with  the  crest  modal  frequencies.  Table  11  summarizes  the  results.  The  cor¬ 
relation  between  most  mode  shapes  was  favorable,  and  the  governing  mode  shape 
in  cross  section  resembled  the  fundamental  mode  shape  of  a  cantilever  beam. 
However,  at  the  lower  frequencies,  the  mode  shapes  for  spillway  monolith  22 
(Figures  52  through  56)  exhibited  discontinuities  at  the  bottom  of  the  pier. 
This  is  probably  because  the  force  input  was  not  large  enough  to  completely 


excite  the  entire  monolith  at  the  lower  frequencies.  As  noted  previously,  the 
spillway  monoliths  are  the  stiffest  part  of  the  dam.  From  a  three-dimensional 
dynamic  FE  analysis  of  an  isolated  spillway  monolith,  the  fundamental  fre¬ 
quency  was  11.3  Hz.  This  is  higher  than  the  first  four  experimentally  mea¬ 
sured  natural  frequencies  of  the  dam  without  reservoir  (5.9  Hz  to  8.2  Hz). 

40.  Comparisons  of  analytical  mo.e  shapes  in  cross  section  with  a 
design  mode  shape  proposed  by  Chopra  (1978)  were  conducted  and  are  presented 
in  Appendix  D.  The  results  were  favorable,  indicating  that  the  proposed 
design  mode  shape  is  sufficiently  accurate  for  design  purposes. 

Dam  With  Reservoir 

41.  Comparisons  of  the  dam  crest  mode  shapes  are  shown  in  Figures  57 
through  61.  Mote  that  the  mode  numbers  correspond  to  bending  modes  of  the  dam 
crest.  Table  12  summarizes  the  results  of  these  comparisons.  (The  response 
at  monolith  22  was  not  used  in  computing  the  MAC  because  excessive  local  re¬ 
sponses  were  recorded  in  test  2  at  this  monolith.)  Crest  mode  shape  correla¬ 
tion  between  experiment  and  analysis  was  not  as  favorable  as  with  the  dam 
without  reservoir.  The  fixed  model  had  higher  MAC  values  than  the  flexible 
base  model  for  the  first  three  modes,  possibly  because  the  prototype  founda¬ 
tion  stiffness  was  greater  than  modeled  in  the  FE  flexible  base  analysis. 

42.  Various  factors  may  have  caused  the  low  correlation  between 
experimental  and  analytical  mode  shapes.  As  discussed  previously,  errors  in 
the  experimental  mode  shapes  may  be  due  partly  to  a  force  input  that  was  not 
large  enough  to  completely  excite  the  entire  dam  at  the  lower  frequencies. 
Also,  a  high  degree  of  modal  coupling  (closely  spaced  resonance  frequencies) 
existed  at  a  higher  frequencies. 

43.  Furthermore,  inadequacies  in  the  FE  model  assumptions  could  have 
produced  a  lower  correlation  between  experimental  and  analytical  mode  shapes: 

a.  The  modified  Westergaard  (1933)  added  mass  concept  to  model  the 
reservoir  effects  may  not  be  valid  at  the  higher  frequencies 
(Bevins,  Chiarito,  and  Hail,  in  preparation).  Westergaard ' s  method 
assumes  that  water  is  incompressible.  This  assumption  is  valid 
when  the  forcing  frequency  is  less  than  the  natural  frequency  of 
the  dam/reservoir  system,  but  for  frequencies  higher  than  this,  the 
method  overest imates  the  hydrodynamic  effects  (Chakrabarti 
and  Chopra  1974).  Comparing  FE  results  with  and  without  reser¬ 
voir,  the  FE  analyses  indicated  that  the  percent  reductions  in 
natural  frequencies  were  greater  for  the  higher  modes.  This  is 
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inconsistent  with  the  experimental  results,  which  indicated  that  j 

the  first  mode  had  the  largest  decrease  in  natural  frequency  due  to  ' 

reservoir  interaction.  t 

i 

b.  The  FE  model  did  not  include  the  effects  of  the  added  mass  of  the 
soil  at  the  abutments.  The  abutment  of  the  Georgia  side  of  the  aam 
was  not  complete  at  the  time  the  first  vibration  test  was  made. 

The  additional  soil  during  the  second  test  might  have  caused  a 
slight  increase  in  the  stiffness  of  the  structure  at  the  Georgia- 
side  abutment. 

c.  Non  1  inear  it  ites  due  to  joint  slippage  were  r.ot  modeled. 

44.  Similar  to  the  case  of  the  dam  without  reservoir,  the  frequencies 
estimated  by  the  FE  models  were  higher  than  measured  experimentally  (7  to 
17  percent  higher  for  the  FE  fixed  base  model,  and  up  to  8  percent  higher  for 
the  FE  flexible  base  model).  The  flexible  base  model  frequencies  compared 
oetter  with  the  experiment  then  did  the  fixed  base  model  frequencies  because 
foundation  interaction  reduces  the  calculated  frequencies. 

45  Comparisons  of  the  experimental  mode  shapes  in  cross  section  with 
results  from  the  FE  model  with  fixed  foundation  are  shown  in  Figures  62 
through  71.  (Comparisons  for  monolith  22  were  not  made  because  only  two  ac¬ 
celerometers  were  mounted  on  monolith  22).  Note  that  the  frequencies  of  the 
cross-sectional  modes  correspond  with  the  crest  modal  frequencies.  Table  13 
summarizes  the  results.  The  correlation  between  most  mode  shapes  was  favor¬ 
able,  and  the  governing  mode  shape  in  cross  section  resembled  the  fundamental 
mode  shape  of  a  cantilever  beam. 

46.  It  should  be  noted  that  some  dam  crest  experimental  mode  shapes 
had  better  correlation  with  FE  results  of  a  different  mode  of  vibration.  The 
fourth  experimental  mode  shape  at  7.5  Hz  resembled  a  fourth  normal  mode  of 
vibration  of  a  beam  fixed  at  both  ends  (Harris  and  Crede  1961).  However,  this 
mode  shape  had  a  better  correlation  with  the  third  FE  crest  flexural  mode. 

Engineering  judgment  was  used  in  interpreting  the  experimental  mode  at  7.5  Hz 
as  a  fourth  mode  of  vibration.  However,  it  is  possible  that  the  experimental 
modes  exhibited  measured  local  behavior  not  determined  by  the  model.  In  this 
way,  the  MAC  can  aid  in  illustrating  complex  responses  not  computed  by  the 
model,  but  measured  in  the  prototype  test. 
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PART  VI :  CONCLUSIONS 


47.  By  visual  inspection  of  the  mode  shapes  with  and  without  reservoir, 
hydrodynamic  interaction  changed  the  mode  shapes  of  the  dam,  with  the  largest 
differences  occurring  at  the  ends  of  the  dam.  The  dam's  midsection  had 
similar  responses  with  and  without  reservoir  interaction;  however,  slight 
shifts  in  the  maximum  response  points  towards  the  ends  of  the  dam  indicate 
that  hydrodynamic  effects  added  mass  to  the  dam's  midsection.  Various  factors 
could  have  caused  fairly  low  correlations  between  experimental  mode  shapes  for 
the  tests  with  and  without  reservoir.  The  force  input  was  not  large  enough  to 
completely  excite  the  entire  structure  at  the  lower  frequencies,  causing 
scatter  in  estimates  of  the  mode  shapes.  Ideally,  two  or  more  force  inputs 
would  be  better.  Also,  a  high  degree  of  modal  coupling  existed  at  the  higher 
frequencies,  and  the  friction  forces  between  monoliths  may  have  been  different 
because  the  tests  were  conducted  during  different  seasons. 

48.  The  experimental  natural  frequencies  were  reduced  by  the  hydro- 
dynamic  effects,  with  the  percent  change  in  the  first  mode  the  largest 

(  10. ‘  percent).  With  the  exception  of  the  fifth  mode,  the  reservoir 
impoundment  had  little  effect  on  the  damping  estimates. 

49.  The  hydrodynamic  pressures  measured  in  the  second  test  with  reser¬ 
voir  were  close  to  the  predicted  values  for  the  fundamental  frequency  (except 
at  elevation  430),  and  the  pressure  distribution  was  similar  to  distributions 
measured  at  Pine  Flat  Dam  (Rea,  Liaw,  and  Chopra  1972).  The  theoretical  hy¬ 
drodynamic  pressures  were  based  on  a  two-dimensional  response  of  an  isolated 
nonoverflow  monolith  (monolith  7).  The  fundamental  mode  shape  of  a  cantiliver 
beam  was  used  to  compute  the  theoretical  pressure  distribution.  Because  the 
experimental  and  analytical  cross-sectional  mode  shapes  of  the  entire  dam 
resembled  the  fundamental  mode  shape  of  a  cantilever  beam,  the  use  of  a  two- 
dimensional  model  was  valid  for  predicting  the  pressure  distribution  on  this 
nonoverflow  monolith.  Further  pressure  measurements  and  analyses  would  be 
needed  to  confirm  if  the  theoretical  pressure  distribution  is  valid  for  the 
other  parts  of  the  dam. 

50.  Measurements  from  both  tests  at  the  relative  joint  motion  arrays 
indicated  relative  motions  between  monoliths.  The  nonlinear  behavoir  of  the 
dam  is  attributed  primarily  to  the  monoliths'  Joint  behavior. 

51.  Despite  a  few  experimental  errors,  no  gross  anomalies  existed,  and 
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Che  dynamic  properties  determined  are  reasonable.  It  can  be  concluded  that 
the  dam  appears  to  be  structurally  sound  as  built. 

52.  Three-dimensional  linear  FE  analyses  have  been  compared  with  previ¬ 
ous  experimental  estimates  of  the  modal  parameters  of  a  concrete  gravity  dam. 
Using  available  dynamic  material  properties  of  the  dam  concrete  and  the  foun¬ 
dation,  a  three-dimensional  FE  model  was  successfully  developed  to  estimate 
the  linear  elastic  modal  properties  of  the  dam.  The  fourth  mode  shape  was  not 
measured  experimentally  for  the  dam  without  reservoir,  but  was  computed  by  the 
three-dimensional  FE  analyses. 

53.  The  MAC  was  useful  for  evaluating  the  correlation  between  the  ex¬ 
perimental  and  the  FE  mode  shapes.  Favorable  correlations  of  the  crest  mode 
shapes  were  computed  for  the  dam  without  reservoir,  indicating  that  a  reason¬ 
able  three-dimensional  FE  model  was  developed  for  computing  the  first  three 
natural  modes  of  vibrations.  Various  factors  may  have  caused  a  lower  cor¬ 
relation  between  the  experimental  and  analytical  crest  mode  shapes  for  the  dam 
with  reservoir.  Low  force  input  at  lower  frequencies  caused  poorer  estimates 
in  experimental  mode  shapes  at  certain  locations,  and  a  high  degree  of  modal 
coupling  existed  at  the  higher  frequencies.  The  effects  of  the  soil-dam  in¬ 
teraction  at  the  abutments  and  the  reservoir  effects  may  not  have  been  modeled 
adequately.  Also,  nonlinearities  due  to  joint  slippage  were  not  modeled.  The 
correlation  between  most  cross-sectional  shapes  was  favorable,  and  the  govern¬ 
ing  mode  shape  in  cross  section  resembled  the  fundamental  mode  shape  of  a 
cantilever  beam. 

54.  It  has  been  shown  that  a  proposed  design  cross-sectional  mode  shape 
by  Chopra  (1978)  is  sufficiently  accurate  for  design  purposes.  Chopra  pre¬ 
sents  a  simplified  analysis  procedure  for  computing  stresses  due  to  earthquake 
ground  motion  in  the  horizontal  direction.  The  simplified  procedure  uses  the 
fundamental  mode  of  vibration  of  a  two-dimensional  model  of  the  dam.  The 
design  mode  shape  resembles  the  fundamental  mode  shape  of  a  cantilever  beam. 

As  discussed  by  Chopra,  the  response  of  concrete  gravity  dams  to  earthquake 
ground  motion  is  primarily  due  to  the  fundamental  mode  of  vibration.  The 
experimental  cross-sectional  mode  shapes  of  the  entire  dam  resembled  Chopra's 
design  mode  shape  (i.e.  the  fundamental  mode  shape  of  a  cantilever  beam). 
Further  studies  would  have  to  be  performed  to  determine  the  validity  of 
Chopra's  method  versus  the  analysis  of  the  entire  dam.  These  further  studies 
would  compare  the  stresses  from  the  analysis  of  the  entire  dam  with  the 
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stresses  from  Chopra's  simpified  procedure. 

55.  Analyses  of  isolated  nonoverflow,  intake,  and  spillway  monoliths 
were  useful  because  they  showed  the  relative  stiffness  of  the  parts  of  the 
dam.  The  calculated  fundamental  frequencies  of  the  isolated  monoliths 
indicated  that  the  intake  sections  were  the  most  flexible  part  of  the  dam  and 
that  the  spillway  monoliths  were  the  stiffest  part  of  the  dam.  Furthermore, 
the  calculated  fundamental  frequency  (7.1  Hz)  of  the  entire  dam  without 
reservoir  for  a  fixed  foundation  was  less  than  the  fundamental  frequencies 

of  the  nonoverflow  and  spillway  monoliths  (9.8  Hz  and  11.3  Hz,  respectively) 
but  slightly  more  than  the  fundamental  frequency  of  the  intake  monolith 
(6.9  Hz).  It  can  be  concluded  from  these  analyses  that  the  intake  sections 
increased  the  flexibility  of  the  entire  dam. 

56.  The  three-dimensional  analyses  estimated  frequencies  up  to  30  per¬ 
cent  higher  than  those  measured  experimentally.  This  is  reasonable  since  the 
FE  models  are  an  approximation  of  the  dam  using  a  discrete  number  of  degrees 
of  freedom.  Due  to  foundation  interaction,  the  FE  mode  with  flexible  foun¬ 
dation  predicted  natural  frequencies  closer  to  the  experimentally  measured 
values  than  did  the  FE  model  with  rigid  foundation. 

57.  In  general,  it  appears  that  three-dimensional  linear  FE  analyses 
can  give  a  good  estimate  of  the  natural  frequencies  of  concrete  gravity  dams 
with  or  without  reservoir  effects.  The  FE  analyses  can  adequately  estimate 
the  bending  mode  shapes  of  concrete  gravity  dams  without  reservoir  impound¬ 
ment.  These  analyses  can  adequately  estimate  the  bending  mode  shapes  of  con¬ 
crete  gravity  dams  with  reservoir  if  the  reservoir  effects  are  appropriately 
modeled . 


PART  VII:  RECOMMENDATIONS  FOR  FURTHER  STUDY 


58.  The  effects  of  the  hydrodynamic  interaction  should  be  studied 
further.  Additional  analyses  could  be  performed  using  a  more  appropriate 
distribution  of  added  masses  due  to  reservoir  impoundment.  The  added  mass 
distribution  should  be  consistent  with  results  from  the  latest  research  on 
hydrodynamic  interaction.  Possibly  a  better  correlation  would  then  be 
achieved  between  analysis  and  the  prototype  test  with  reservoir  impoundment. 

59.  Other  refinements  in  the  FE  analysis  should  include  more  accurate 
modeling  of  soil-structure  interaction  and  a  model  for  the  joint  slippage 
between  monoliths. 

60.  Studies  should  be  performed  on  the  data  collected  measuring  the 
movement  of  the  dam  before  and  after  reservoir  impoundment.  This  would  give 
an  estimate  of  the  flexibility  of  the  interacting  dam-foundation  system. 

61.  Further  pressure  measurements  and  analyses  would  be  needed  to 
confirm  if  the  theoretical  pressure  distribution  (Fenves  and  Chopra  1986)  is 
valid  for  the  entire  dam.  Also  further  studies  would  have  to  be  performed  to 
determine  the  validity  of  Chopra's  simplified  method  of  analysis  (Chopra  1978) 
versus  the  analysis  of  the  entire  dam.  These  further  studies  would  compare 
the  stresses  from  the  analysis  of  the  entire  dam  with  the  stresses  from 
Chopra's  simplified  procedure. 
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Table  1 

Results  from  Uniaxial  Compression  Test 


Cy 1 inder 
No. 


Poisson ' s 
Ratio 


Initial 

Elastic 

Concrete 

Modulus 

1 n6  ; 


Uniaxial 

Compressive 

Strength 

ksi 


Date* 

Tested 


5,037 
5,81 1 


4  ksi  Design  Strength  Concrete 


3  ksi  Design  Strength  Concrete 


=  6  May  1982,  2  =  27  Aug  1989. 


5  Mar 
10  Jul 


5 , 464 

0.18 

4.51 

5.68 

8  May 

5,765 

0.16 

3-45 

5.20 

2  Jul 

t> ,  468 

0.16 

3.76 

4.55 

9  May 

5,563 

0.18 

3.96 

4.95 

29  May 

5,735 

0. 16 

4.09 

4.85 

26  Jun 

6,112 

0.16 

4.90 

6.16 

19  Sep 

6,140 

0.19 

4.45 

5.22 

25  Sep 

5,166 

0.20 

4.91 

5.86 

20  Mar 

5,512 

0.16 

3.70 

4.81 

14  May 

5,655 

0.13 

3-92 

5.25 

12  Jun 

5,708 

0.15 

4.12 

4.81 

19  Jun 

5,214 

0. 16 

4.02 

4.72 

25  Mar 

5,519 

0.18 

4.46 

5.11 

15  May 

6,016 

0.15 

4.22 

6.30 

27  Aug 

6,083 

0.18 

4.92 

7.27 

11  Sep 

6,187 

0.13 

3.24 

3.87 

1  Oct 

2  ksi  Design  Strenf 

Sth  Concrete 

4,954 

0.20 

3.44 

3.45 

25  Feb 

5,833 

0.  19 

3.50 

3.40 

17  Jul 

4,993 

0. 18 

3.45 

3.13 

28  Feb 

5,993 

0. 18 

2.63 

3.  17 

21  Aug 

m 


85» 


Table  2 

Initial  Elastic  Concrete  Modulus 


Cylinder 

No. 


5,037 

5,811 


^static 

c 

F 

dynamic 

c 

Percent 

Date* 

x  10°  psi* 

x  10°  psi* 

Increase 

Tested 

4  ksi  Design  Strength  Concrete 


4.67 

5.59 

20 

1 

4.12 

4.61 

12 

2 

16  Average 


3  ksi  Design  Strength  Concrete 


5,765 

3.45 

4.72 

72 

5,468 

3.76 

4.48 

19 

5,563 

3-96 

4.63 

17 

5,735 

4.09 

4.72 

15 

6,112 

4.90 

5.64 

15 

6,140 

4.45 

4.92 

1 1 

5,512 

3.70 

5.34 

44 

5,655 

3.92 

4.81 

23 

5,708 

4.12 

4.77 

16 

5,214 

4.02 

5.08 

26 

5,519 

4.46 

5.08 

14 

6 , 0 1 6 

4.22 

5.48 

30 

6,083 

4.92 

5.73 

16 

6,187 

3.24 

4.92 

52 

5,464 

4.51 

5.43 

20 

5,166 

4.91 

5.81 

18 

23  Average 


2_ 

ksi  Design  Strength  Concrete 

4,993 

3.45 

3.85 

12 

5,993 

2.63 

4.09 

56 

4,954 

3.44 

4.54 

32 

5,833 

3.50 

4.54 

29 

2 


1 

1 


2 

2 

1 

1 


32  Average 
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Table  4 

Dynamic  Modulus  of  Elasticity  of 
Rock  Core  Samples 


o] 


Modulus  of  Elasticity 


Static 
( Initial 

Tangent) 

Dynamic* 

7.14 

0.54f 

5.20 

5.78 

12.90 

0.40t 

5.20 

7.43 

13.79 

1 .05t 

3.18 

3.90 

5.97 

0.24f 

7.63 

9.60 

4.53 

5.70 

7  sample  average 
3  sample  average 


The  dynamic  modulus  was  estimated  from  the  three-sample  average  by  apply 
ing  the  same  ratio  of  dynamic  to  static  modulus. 

9.60  =  x  7.63 

Not  used  in  average.  Sample  contained  healed  fractures  that  probably 
interfered  with  frequency  reading. 


Table  5 

Experimental  Damping  Estimates  Before  and  After 


Reservoir  Impoundment 


Damping 

Damping 

Before 

After 

Impoundment 

Impoundment 

percent 

percent 

Not  measured 


Table  6 

Summary  of  Experimental  and  Analytical  Frequencies,  Hz 


Crest 

Mode 

No. 

1 

2 

3 

4 

5 


Crest 
Mode 
No . 


Prototype 
Vibration  Test 

Finite  Element 
Analysis  Without 
Reservoir 

Without 

Reservoir 

Fixed 

Base 

Flexible 

Base 

5.9 

7.1 

6.6 

6.8 

7.9 

7.4 

7.6 

9.3 

8.7 

Not  measured 

10.4 

9.6 

8.2 

10.7 

10.0 

Finite  Element 


Prototype 

Analysis  With 

Vibration  Test 

Reservoir 

With 

Fixed 

Flexible 

Reservoir 

Base 

Base 

5.3 

6.2 

5.7 

6.2 

6.8 

6.4 

7.3 

7.8 

7.3 

7.5 

8.2 

7.6 

8.1 

8.7 

8.0 

Table  7 

Experimental  Mode  Shape  Comparisons 


Prototype 
Vibration  Test 
Without 
Reservoir 
Frequency  Hz 


Protype 

Vibration  Test 
With 

Reservoir 
Frequency,  Hz 


MAC* 


Table  8 

Finite  Element  Fixed  Base  Mode  Shape  Comparisons 


Crest 

Mode 

No. 


Crest 

Mode 

No. 


Analysis 
Without 
Reservoir 
Frequency,  Hz 


analysis 
Without 
Reservoir 
Frequency,  Hz 


Analysis 

With 

Reservoir 
Frequency,  Hz 


Analysis 

With 

Reservoir 
Frequency.  Hz 


Table  10 

Crest  Mode  Comparisons  -  Dam  Without  Reservoir 


(Experiment  versus  Analysis) 


Crest 

Mode 

No. 


Prototype 

Vibration 

Test 


Three-Dimensional  Finite  Element  Analyses 


F  ixed 
Base 


Flexible 

Base 


Frequency,  Hz 

Frequency,  Hz 

MAC* 

Frequency,  Hz 

MAC* 

5.9 

7.1 

0.73 

6.6 

0.74 

6.8 

7.9 

0.63 

7.4 

0.63 

7.6 

9-3 

0.64 

8.7 

0.61 

Not  measured 

10.4 

— 

9.6 

-- 

8.2 

10.7 

0.41 

10.0 

0.  13 

Table  11 

Cross  Section  Comparisons  Without  Reservoir 
(Experiment  versus  Analysis) 


Crest 
Mode 
No . 


Experimental 
Frequency,  Hz 

5.9 

6.8 

7.6 

Not  measured 


Finite  Element 
Fixed  Base 
Frequency,  Hz 


Frequenc\ 

7.1 

7.9 

9.3 

10.4 


Modal 

Monolith 

7 


Assurance 
Monoi ith 
16 


Criterion 

Monolith 

22 
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Figure  3.  Crest  accelerometer  locations  for  test  1 
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Figure  4.  Crest  accelerometer  locations  for  test  2 
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Figure  6.  Cross-sectional  accelerometer  and  pressure  gage  locations  for  test  2 
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Figure  7.  Relative  Joint  motion  accelerometer  locations 
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Figure  11.  Experimental  mode  shape  4  comparisons 
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Figure  14.  Comparisons  of  experimentally  measured  and 
theoretical  hydrodynamic  pressures 
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Figure  19.  Eigenvectors  Nos.  6,  8,  9,  11,  and  14  from  FE  model 
without  reservoir  and  with  fixed  base 


EIGENVECTOR  NO.  1.  6.2  Hz  (CREST  MODE  NO.  1) 


i  £ .  a'  a.  .--  .-•r 


EIGENVECTOR  NO.  2,  6.3  Hz 


/  /ri .  5  ; 


EIGENVECTOR  NO.  3,  6.6  Hz 


EIGENVECTOR  NO.  11  6.8  Hz  (CREST  MODE  NO.  2) 


EIGENVECTOR  NO.  14.  7.8  Hz  (CREST  MODE  NO.  3) 

Figure  20.  Eigenvectors  Nos.  1,  2,  3,  11,  and  14  from  FE  model 
with  reservoir  and  fixed  base 
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Figure  21.  Eigenvectors  Nos.  16-20  from  FE  model  with 
reservoir  and  fixed  base 
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Figure  22.  FE  fixed  base  model  mode  shape  1  comparisons 
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Figure  27.  FE  flexible  base  model  mode  shape  1  comparisons 
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Figure  39.  Cross-sectional  mode  shape  3  comparisons  For 
monolith  7,  upstream  face,  dam  without  reservoir 
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Figure  41.  Cross-sectional  mode  shape  5  comparisons  for 
monolith  7,  upstream  face,  dam  without  reservoir 
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Figure  48.  Cross-sectional  mode  shape  2  comparisons  For 
monolith  1 6 ,  dam  without  reservoir 
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monolith  22,  dam  without  reservoir 
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Figure  55.  Cross-sectional  mode  shape  4  comparisons  for 
monolith  22,  dam  without  reservoir 
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Figure  66.  Cross-sectional  mode  shape  5  comparisons  for 
monolith  7,  downstream  face,  dam  with  reservoir 
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Figure  70.  Cross-sectional  mode  shape  4  comparisons  for 
monolith  16,  dam  with  reservoir 
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Figure  71.  Cross-sectional  mode  shape  5  comparisons  for 
monolith  1 6 ,  downstream  face,  dam  with  reservoir 


APPENDIX  A:  FOUNDATION  INTERACTION  STUDIES 


1.  Foundation  properties  are  taken  from  Richard  B.  Russell  Dam  Design 
Memorandum  8,  1983.  The  in  situ  elastic  modulus  of  the  foundation,  E^  ,  was 
estimated  from  the  average  dynamic  modulus  of  elasticity,  cross-hole  seismic 
velocity  and  pulse  velocity  of  rock  core  samples.  The  dynamic  modulus  of 
elasticity,  ERq  ,  measured  from  seven  rock  core  samples  is  shown  in  Table  4 
of  the  main  text. 


Erc  =  9.6  x  106  psi 

2.  The  average  cross-hole  seismic  velocity,  v^.  ,  taken  from  44  mea¬ 
surements  was  12,810  ft/sec,  with  a  coefficient  of  variation,  Cy  =  0.22  . 
The  average  pulse  velocity,  VR£  ,  taken  from  10  rock  core  samples  was 
16,130  ft/sec  with  Cv  =  0.20  .  The  velocity  of  waves,  v  ,  propagating 
through  an  elastic  material  is  (Timoshenko  and  Goodier  1970)* 

where 

E  =  modulus  of  elasticity 
0  =  mass  density 

Assuming  constant  mass  density,  Ef 
tion  A1 : 


E  "2 
RC  VRC 

Then 

Ef  =  0.63  x  9.6  x  106  psi 
=  6.05  x  106  psi 

3-  Elastic  theory  was  used  to  obtain  spring  constants  for  the  flexible 
base  analysis  of  the  entire  dam.  The  stiffness  of  a  semi- i nf in i te  solid  acted 
on  by  a  distributed  loading  on  the  boundary  is  (Timoshenko  and  Goodier  1970) 


V 


(A1) 


can  be  estimated  from  ERC  using  Equa- 


=  0.63 

16,130^ 


(A2) 


*  S^e  .'eferences  at  the  end  of  the  main  text. 

A 1 


■r. 


where 


k  =  stiffness 

E  =  modulus  of  elasticity 

A  =  magnitude  of  distributed  load  area 

m  =  0.95  for  a  square  loading 

=  0.92  for  a  2:1  rectangular  loading 

v  =  Poisson's  ratio 

Using  m  =  0.92  ,  Poisson's  ratio  =  0.2  ,  and  an  elastic  modulus  of  foundation 
of  6.05  x  10®  psi,  spring  constants  were  calculated  for  each  node  along  the 
base  of  the  dam.  Because  the  loaded  area  varied,  the  constants  varied  from 
283  x  10®  lb/ft  to  461  x  10®  lb/ft.  These  constants  gave  an  upper  bound  to 
results  obtained  from  vibration  tests. 

4.  Studies  were  performed  on  the  effect  of  foundation  interaction 
on  the  natural  frequencies  of  individual  two-dimensional  ( 2— D )  and  three- 
dimensional  (3-D)  nonoverflow  sections.  Sections  used  for  the  2-D  and  3-D  FE 
analyses  were  1  and  48  ft  wide,  respectively.  Tables  A1  and  A2  summarize  com¬ 
parisons  of  natural  frequencies  of  the  nonoverflow  section  using: 

a.  Winkler  springs  along  base. 

b.  FE  foundation  (Figure  A1). 

c.  Fixed  foundation. 

5.  The  grid  used  for  the  nonoverflow  section  in  the  analysis  of  the 
entire  dam  was  also  used  for  the  individual  nonoverflow  section  analyses  with 
Winkler  springs  and  fixed  foundation.  The  grid  shown  in  Figure  A1  was  used 
for  the  analysis  with  FE  foundation.  Dynamic  concrete  modulus  <">f  elasticity 
values  were  used,  along  with  an  elastic  modulus  of  foundation  equal  to  9.60 

x  10®  psi.  Figures  A 2  and  A3  show  the  effect  of  spring  stiffness  on  the 
fundamental  frequency  of  the  nonoverflow  monolith. 

6.  The  results  indicated  that  a  relatively  low  spring  constant  approxi¬ 
mated  the  FE  analysis  with  foundation  (Figure  A 2).  For  the  2-D  model,  the  ef¬ 
fect  of  the  horizontal  spring  constant  on  natural  frequency  is  insignificant 

O 

when  it  exceeds  about  2  x  10°  lb/ft,  as  shown  in  Table  A1  and  Figure  A3- 

7.  As  shown  in  Figure  A2,  the  fundamental  frequency  of  the  2-D  model 
approaches  the  fundamental  frequency  with  a  fixed  foundation  for  values  of 

O 

spring  constant  above  20  *  10°  Lb/ft.  However,  the  fundament.il  frequency  of 


fc-wx'vi  IV'  W  '"V  rv  W*  ^  *  ’’’ 


the  3-D  model  approaches  the  fundamental  frequency  with  a  fixed  foundation  for 

Q 

values  of  spring  constant  over  100  x  10  lb/ft.  This  indicates  that  the  3-D 
model  is  more  flexible  than  the  2-D  model. 

8.  Elastic  theory  predicts  a  much  stiffer  spring  constant  than  obtained 
from  comparison  with  thu  results  using  an  FE  foundation  grid.  Using  a  modulus 
of  elasticity  of  9-60  x  10®  psi,  Equation  A3  gave  a  spring  constant  of  92.7 
x  10®  lb/ft  for  the  2-D  model  and  a  spring  constant  of  383  *  10®  lb/ft  for  the 
3-D  model.  The  equivalent  spring  constants  of  the  FE  foundation  grid  results 
were  4.5  x  10®  lb/ft  for  the  2-D  model  and  12.3  *  10®  lb/ft  for  the  3-D  model 
(Figure  A2).  Thus,  elastic  theory  predicted  foundation  stiffness  20  to  30 
times  larger  than  an  FE  model  with  foundation  grid. 
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lement  nonoverflow  monol 


Effect  of  horizontal  spring  stiffness  on  fundamental  frequency  of 
nonoverflow  monolith  (1  lb/ft  =  14.6  N/m) 


1.  The  modal  assurance  criterion  (MAC)  is  a  measure  of  the  consistency 
between  estimates  of  modal  vector.  This  criterion  was  used  to  indicate  quan¬ 
titatively  how  well  the  mode  shapes  compared.  The  modal  assurance  criterion 
is  defined  as 


MAC(c  d)  -  - 1  M0M( c , d )  | - 

mauc.cj;  -  M0M(C)C)  M0M(d ,d) 


where 

MAC(c,d)  =  MAC  for  vectors  c  and  d 

M0M(c,d)  =  cross  moment  of  the  modal  vectors  c  and  d 
MOM(c,c)  =  auto  moment  of  the  mode  vector  c 
M0M(d,d)  =  auto  moment  of  the  mode  vector  d 
For  real  modes.  Equation  B1  reduces  to 

MAC 

where 


( B2 ) 


$.  -  element  of  mode  vector  <& 

0.  =  element  of  mode  vector  0 

i 

MAC  =  MAC  of  vectors  $  and  0 

2.  The  MAC  criterion  will  have  values  from  zero,  representing  no  cor¬ 
respondence,  to  one,  representing  a  consistent  correspondence. 


*  Allemang  and  Brown  (1982).  See  References  at  the  end  of  the  main  text. 
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APPENDIX  C:  HYDRODYNAMIC  INTERACTION 

1.  Westergaard  (1933)*  derived  an  approximate  formula  for  the  added 
mass  of  water  on  the  vertical  upstream  face  of  a  dam  during  earthquakes.  As 
shown  in  Figure  Cl,  Westergaard  assumed  a  parabolic  added  mass  distribution. 
Using  consistent  units: 


b 


(Cl) 


where 

b  =  width  of  added  water  mass 
h  r  depth  of  reservoir 
y  =  vertical  distance  from  top  of  dam 
Westergaard  neglected  the  effects  of  flexibility  of  the  dam  and  compress¬ 
ibility  of  water. 

2.  Chopra  (1978)  developed  an  approximate  analysis  method  accounting 
for  the  flexibility  of  the  dam  and  compressibility  of  water.  The  approximate 
natural  period  of  vibration  of  the  dam  with  reservoir  effects  is 


T 


s 


R1  1.4 


( C2 ) 


where 

Ts  =  vibration  period  with  water,  sec 

R i  =  ratio  of  fundamental  vibration  periods  of  dam  with  and 
without  water,  plotted  (Chopra  1978)  against  depth  of 
water  for  various  values  of  concrete  modulus  of  elasticity 

Hs  =  height  of  dam,  ft 
E  -  concrete  modulus  of  elasticity,  psi 

3.  For  the  Richard  B.  Russell  Dam,  assume  that  height  of  dam  -  185  ft, 
depth  of  water  =  170  ft,  and  concrete  modulus  of  elasticity  =  5  *  10^  psi. 
Then  P1  =  1.31  ,  and  the  vibration  period  computed  by  Equation  C2  is 
0.’517  sec.  Therefore,  the  fundamental  natural  frequency  of  the  dam  with 
water  is  6.59  Hz  using  Chopra's  approximate  analysis. 


*  See  References  at  the  end  of  the  main  text. 


4.  As  discussed  below,  a  modified  Westergaard  formula  was  used  to  com¬ 
pute  a  fundamental  natural  frequency  nearly  the  same  as  that  computed  using 
Chopra's  procedure.  As  shown  in  Figure  Cl,  the  mass  distribution  from 
Equation  Cl  was  factored  by  0.5,  producing  a  total  water  force  of  43,832  lb  on 
a  1-in. -wide  face  of  dam: 


b  = 


5.  As  shown  in  Figure  Cl,  the  constant  reservoir  width,  W  ,  that 
produced  the  total  water  force  from  Equation  C3  was  calculated. 

total  water  force 


where 

W  =  constant  reservoir  width 

Yw  =  weight  density  of  water 
h  =  depth  of  reservoir 

With  y  =  62.4  lb/ft^  ,  h  =  170  ft  ,  and  total  water  force  =  43,832  lb/in., 
w 

the  constant  reservoir  width,  W  ,  from  Equation  C4  is  49.6  ft.  A  fundamental 
natural  frequency  of  6.62  Hz  was  obtained  from  a  two-dimensional  dynamic 
finite  element  (FE)  analysis  of  the  section  shown  in  Figure  Cl  with  concen¬ 
trated  masses  applied  to  the  nodes  in  contact  with  the  constant  49.6-ft  reser¬ 
voir  width.  The  same  dimensions  and  material  properties  of  the  section  used 
in  Chopra's  analysis  were  used  in  this  FE  analysis. 

6.  The  fundamental  natural  frequency  of  dam  with  reservoir  calculated 
from  the  uniform  width  of  reservoir  (6.62  Hz)  compared  well  with  the  result 
from  Chopra  (1978)  (6.59  Hz,  or  0.5-percent  difference).  Therefore,  the  uni¬ 
form  reservoir  width  concept  was  used  to  obtain  added  masses  on  the  model  of 
the  entire  dam.  As  a  first  approximation,  a  constant  49.6-ft-wide  mass  of 
reservoir  was  applied  at  all  sections  of  the  dam.  Concentrated  masses  due  to 
the  reservoir  were  applied  at  nodes  assumed  to  be  in  contact  with  the 
reservoir . 

7.  Additional  studies  on  hydrodynamic  interaction  would  use  Chopra's 
mass  distribution  in  the  vertical  direction  and  a  modified  nonuniform  mass 
distribution  in  horizontal  direction  along  the  length  of  the  dam. 
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APPENDIX  D:  COMPARISON  OF  ANALYTICAL  AND  DESIGN 
CROSS-SECTIONAL  MODES 


1.  The  modal  assurance  criterion  (MAC)  (Allemang  and  Brown  1982)*  was 
used  to  quantitatively  compare  mode  shapes  in  cross  section  from  finite 
element  (FE)  results  with  the  design  mode  shape  proposed  by  Chopra  (1978). 

2.  In  Figures  D1  to  D3,  the  design  mode  shape  is  compared  with  three 
mode  shapes  from  three-dimensional  (3-D)  FE  analyses  to  the  Richard  B.  Russell 
Dam  without  reservoir,  fixed  base,  for  nonoverflow  monolith  6,  spillway  mono¬ 
lith  *7,  and  intake  monolith  11,  respectively.  The  three  FE  modes  shapes  are 
at  the  normal  frequencies  corresponding  to  the  first  three  FE  normal  mode 
shapes  of  the  dam  crest.  All  shapes  resemble  the  fundamental  mode  shape  of  a 
cantilever  beam.  As  summarized  in  Table  D1,  the  high  value  of  the  MAC  in¬ 
dicates  that  the  analytical  and  design  mode  shapes  compare  well. 

3.  The  analytical  mode  shapes  compared  well  with  the  experimental  mode 
shapes  in  cross  section  as  stated  in  the  report.  It  can  be  concluded  that  the 
design  mode  shape  would  also  compare  well  with  the  experimental  results. 
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Figure  D2.  Analytical  and  design  mode  shape 
comparisons  For  spillway  monolith  17 
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Figure  D3.  Analytical  and  design  mode  shape 
comparisons  for  intake  monolith  11 


APPENDIX  C:  GLOSSARY 


Acce 1  er  at  ior.  -  A  vector  quantity  that  indicates  the  time  rate  of  change  of 

ve oe  i  ty  . 

Coherence  Foret  ion  -  The  coherence  function  for  an  input/output  pair  measures 
what  ratio  of  the  output  acceleration  is  I  inear ly  caused  by  a  measured  input 

force . 

Coup iec  Modes  -  Modes  of  vibration  that  are  dependent  upon  other  modes  of 
vibration  due  to  energy  transfer.  (See  "Mode  of  vibration.") 

Critical  Damping  -  The  minimum  amount  of  damping  for  which  no  oscillation 
occurs  when  a  displaced  system  returns  to  its  original  position. 

Cycle  -  A  set  of  values  that  occurs  during  one  complete  performance  of  a 
periodic  process. 

Damped  Natural  Frequency  -  Frequency  of  free  vibration  of  a  damped  linear 
system. 

Damp  ing  -  Energy  loss  of  a  system  due  to  friction  and  other  resistances. 

Damping  Ratio  -  Ratio  of  the  actual  amount  of  damping  of  a  system  to  the 
critical  damping  of  the  system. 

Degrees-of -freedom  -  The  number  of  independent  coordinates  required  to  define 
the  motion  of  a  system. 

Dynamic  Modulus  of  Elasticity  -  A  measure  of  the  material  elasticity  under 
dynamic  loads,  calculated  from  the  measurement  of  the  fundamental  resonant 
frequency  of  vibration. 

E igenval ue  -  Frequency  of  a  normal  mode  of  vibration. 

E i genvector  -  A  configuration  of  a  vibrating  system,  or  normal  mode  of 
vibration  in  which  ail  particles  are  in  harmonic  or  sinusoidal  motion  at  the 
same  frequency.  Each  eigenvector  can  exist  independently  of  other  eigen¬ 
vectors  of  the  system. 

Excitation  -  An  external  force  acting  on  a  system  that  causes  the  system  to 
rest'  id  in  a  particular  way. 

Forces  Vibration  -  Vibration  or  oscillation  of  a  system  due  to  excitation 
( external  forces ) . 


Free  Vioration  -  Oscillation  of  a  system  under  the  action  of  internal  forces 
of  the  system,  without  external  excitation. 

Fr  ■  [  a  q  .-y  -  The  roc.’  procal  of  the  period  or  repetition  time  of  oscillation, 
ex:u  cr-sod  in  eye.-*  per  unit  time.  Hertz  (Hz)  is  cycle  per  second. 

- r •  'T-  uonse  Fur.ot  i  )t;  :  Tr  ansfer  Function)  -  The  complex  (  i . e . ,  composed 

;  f  ■  a  i  sn  i  ;  mag  it.  try  parts)  ratio  of  the  measured  output  acceleration  to  the 
r o  I  ’  ue  . 


it.-'*!,  -roe  .•••icy  -  The  lowest  natural  frequency  of  a  system. 

iar.-r.t  i .  Mode  of  7:  nr  tt,  ion  -  Mode  that,  has  the  lowest  natural  frequency. 
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